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A  semi-empirical  method  involving  asymptotic  expansions  is  used  to  obtain  an  approxi¬ 
mate  formula  for  the  fundamental  frequency  of  a  uniform  rotating  beam  clamped  off  the 
axis  of  rotation.  Results  from  the  formula  are  shown  to  be  of  the  order  of  0-1%  different 
from  the  exact  results  for  a  wide  range  of  rotor  speeds  and  hub  radii  up  to  the  order  of  blade 
length.  Thus,  the  designer  is  provided  with  a  rapid,  very  accurate  estimate  of  the  frequency, 
without  having  to  interpolate  results  from  a  chart  or  run  a  digital  computer  program. 

1.  INTRODUCTION 

The  problem  of  approximating  the  fundamental  frequency  of  a  uniform  rotating  beam  has 
application  in  the  design  of  helicopter  rotor  blades,  wind  turbine  blades,  and  flexible 
satellite  booms.  Although  the  Ritz  method  [1]  can  yield  quite  accurate  estimates  of  the 
fundamental  frequendy^fvith  a  moderate  number  of  terms,  it  is  sometimes  desirable  to 
have  a  fairly  accurate  formula  that  can  be  evaluated  by  using  a  pocket  calculator  rather 
than  having  to  interpolate  from  charts  or  run  a  digital  computer  program.  Such  an 
expression  is  not  trivial  to  come  by,  however,  since  the  exact  solution  to  the  problem  exists 
only  as  an  infinite  series  [2].  Uniform  approximations  based  on  a  small  parameter,  such  as 
17  =  EI/mfi2L 4,  are  accurate  only  when  file  parameter  is  quite  small  and  are  thus  not 
general  enough  for  practical  use  [3].  A  useful  approach  appears  to  be  a  semi-empirical 
method  referred  to  as  a  “composite”  expansion,  as  developed  by  Peters  [4]  for  rotating 
beams  clamped  at  the  axis  of  rotation.  In  Peters'  work  [4]  the  expansion  was  accurate  for 
all  values  of  rj,  but  did  not  include  variation  in  »ff -clamping  a  =  R/L  with  R  being  the  hub 
radius.  It  is  the  purpose  of  this  paper  to  extend  the  results  of  reference  [4]  for  the 
fundamental  bending  frequency  to  include  variation  in  the  off-clamping  parameter. 

To  proceed  with  this  expansion,  one  first  notes  that  the  ordinary  differential  equation  for 
the  mode  shape  of  a  rotating  beam  when  the  stiffness  parameter  rj  vanishes  corresponds 
physically  to  that  of  a  rotating  string  with  its  root  offset  from  the  axis  of  rotation.  An 
approximate  composite  expansion  solution  for  the  fundamental  frequency  of  the  rotating 
string  can  then  be  used  to  extend  the  composite  expansion  of  reference  [4]  to  include  the 
effect  of  off-clamping. 

2.  STATEMENT  OF  THE  PROBLEM 

The  differential  equation  for  free  vibration?  in  bending  of  a  uniform  rotating  beam  in  a 
plane  at  angle  9  with  the  rotation  axis  (safftngure  1)  is 

tju""-[J(l  -x)u']'~iiu  =0,  (1) 
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Plane  in  which  u 
deflection  occurs 


( a )  (b  ) 

Figure  1.  Schematic  of  rotating  beam  clamped  off  the  axis  of  rotation. 


where 


tj  =EI/mfi2Li,  a=R/L ,  u  =  u/L ,  x  =  x/L,  Os*  si, 

2,~2  .  2  J 0  =  0  out-of -plane bending] 

m=w/722  +  sin20  ,  .  (2) 

1 0  =  7t/2  in-plane  bending  ) 

The  differential  equation  is  an  eigenvalue  problem  in  n  for  which  approximations  are 
sought  as  t)  and  a  are  varied.  The  boundary  conditions  for  a  cantilevered  root  end  are 

i/(0)  =  0,  t/'(0)  =  0,  «"(1)  =  0,  «'"(1)  =  0.  (3) 

One  can  expand  the  quantities  u  and  ji  in  power  series  of  the  small  quantity  tj1/2  as  in 
reference  [4]  so  that  u  =  u0  +  uxrj 1/2  +  ■  •  •  and  ji  =  jix)  +  hxtjU2  +  ■  •  The  equation 
corresponding  to  the  zeroth  power  of  tj  is 

Mo«o  +  {[«U -x)  +  i(\ -X2)]ni}'  =  0.  (4) 

For  a  =  0  it  was  shown  in  reference  [4]  that  only  the  displacement  boundary  condition 
affects  /x  to  order  tj.  Thus, 

wo(0)  =  0.  (5) 

It  is  assumed  that  equation  (5)  holds  for  the  case  when  a^O  as  well.  Equation  (4) 
represents  the  differential  equation  for  a  rotating  string  with  the  root-end  offset  parameter 
a.  It  is  necessary  to  obtain  a  general  approximation  for  n»  (i.e.,  a  small  rj  approximation  for 
n )  from  this  equation  which  will  serve  as  one  ingredient  in  the  process  of  approximating  jx 
for  equation  (1).  Note  that  an  approximation  for  ji  with  17  small  and  large  compared  to 
unity  with  a  =  0  has  already  been  obtained  in  reference  [4].  Here  the  case  for  large  tj  and 
a  /  0  will  be  examined  in  a  later  section. 

The  exact  solution  of  equation  (4)  can  be  expressed  in  terms  of  hypergeometric 
functions.  A  result  more  useful  for  present  purposes  can  be  obtained,  however,  by  looking 
at  limiting  behavior  for  small  and  large  a,  respectively.  To  this  end,  it  is  convenient  to 
rewrite  equation  (4)  temporarily  with  A  =  and  v  =  m0: 

Ac+{[«(1-x)  +  |(1-x2)]c'}'  =  0,  i>(0)  =  0.  (6) 
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3.  ROTATING  STRING  WITH  SMALL  OFFSET 

For  small  a  the  variables  A  and  v  can  be  expanded  in  powers  of  a,  A  =  A()  +  A  i a  +  ■  •  ■ 
and  v  =  c()+  Eitr  +  ■  •  • ,  which  yields  the  following  equations  that  correspond  to  powers  of 
a  from  equation  (6): 

a":  AuUo  +  Ud  -xJ)i’o]'  =  0,  (7a) 

a1:  A«Ui  +![(1  - x2)v\  ]'  =  -A,t>0-[(1  -x)e«]'.  (7b) 

From  reference  [4]  the  solution  is,  for  c<>(0)  =  0,  the  odd  Legendre  polynomials.  For  the 
first  mode  (n  =  1) 

t’o=P2n-i(T)  =  Jt,  A0  =  n(2n  -  1)  =  1.  (8) 

To  obtain  A  i,  it  is  observed  that  v„  must  be  orthogonal  to  the  right-hand  side  of  equation 
(7b).  Thus, 

f  {A,u0  +  [(1 -x)U()]'}t)odx  =  0,  (9) 

Jo 

or 

A i  =  (  (1  -x)v'o  dx/f  Dodx  =  i  (10) 

Jo  /  Jo 

Since  an  estimate  of  the  fundamental  frequency  is  all  that  is  of  interest  here,  one  can  take 

A  =  1  +2«  (ID 

as  an  indication  of  the  behavior  of  A  for  equation  (6)  with  a  small.  It  should  be  noted  that 
the  linear  term  A  i  =  i  does  not  agree  with  A  i  =  6/ir  obtained  in  reference  [3].  The  reasons 
for  this  disagreement  are  not  known;  however,  results  from  an  exact  numerical  calculation 
(see  below)  seem  to  indicate  that  A !  =  |  is  correct. 


4.  ROTATING  STRING  WITH  LARGE  OFFSET 
The  solution  of  equation  (6)  for  large  a  is  similarly  obtained  with  the  definition  A/a  =  v 


and  f  =  1/a  so  that  equation  (6)  becomes 

w  +  {[( 1  -  x )  +  (f/2)(  1  -  *2)]r'}’  =  0,  (12) 

and  one  assumes 

v  —  Vo  +  Ft'i ,  v  =  vn  +  ev\,  t>o(0)  =  ui(0)  =  0.  (13) 

The  equations  for  zeroth  and  first  powers  of  f  are 

f":  v oV„  +  [(1  -*)!>,',]'  =  0,  (14a) 

f1:  FoUi+[(l -mc1]'=  (14b) 

The  solution  for  v„  involves  a  zeroth  order  Bessel  function 

t-„  =  J„[2vVo(l-7)]=  I  (~1)" -/I  — r  (15) 

n-D  (n!) 


The  boundary  condition  rn(0)  =  0  yields 


„„=  1-445796491. 


(16) 
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To  avoid  secular  terms,  v0  must  be  orthogonal  to  the  right-hand  side  of  equation  (14b)  so 
that 

v,  =jj  (l-jt2)i>oZd*/!  no  dx  =  1  038163742.  (17) 


Thus  for  large  a 

A  =  av  =  a(i>o  +  vxe)  =  «i'o+  v\,  (18) 

with  v0  and  v\  given  in  equations  (16)  and  (17),  respectively.  It  is  interesting,  upon 
comparison  of  equations  (11)  and  (18),  to  note  that  the  behavior  for  large  and  small  a  is 
quite  similar.  Therefore,  a  simple  exponential  patching  to  produce  an  accurate  approxi¬ 
mation  of  A  for  all  a  will  now  be  constructed. 


5.  COMPOSITE  EXPANSION  FOR  THE  FUNDAMENTAL  FREQUENCY 
OF  A  ROTATING  STRING 

If  it  is  assumed  that  the  behavior  of  A  for  all  a  can  be  expressed  in  the  form 


A  =  av0+  vi+e  "( a  +  ba ), 

(19) 

then,  for  small  a. 

A  =  av  0+  vx  +  (1  —  a  +  •  •  •)(“  +  ba). 

(20) 

For  small  a,  A  =  A0 

+  A,<*  so  that 

a-  Ao~  pi.  b  ~  A,  -  i/0~  v\  +A0. 

(21) 

Thus,  for  small  tj,  the  smallest  eigenvalue  of  equation  (4)  is  approximately 

Ho  =  avo+vi  +  g(a), 

(22) 

where 

g(a  )  =  e-a[Ao  —  v,  +  (A0  +  A  i  -  v0  -  vt )a]. 

(23) 

Table  1 


Ho  versus  a 


a 

Ha, 

expansion  from 
equation  (22) 

Ha, 

exact 

0 

1 

1 

1 

2-47582 

2-46802 

2 

3-92893 

3-91965 

3 

5-37605 

5-36815 

4 

6-82183 

6-81551 

5 

8-26743 

8-26232 

6 

9-71309 

9-70883 

7 

11-1588 

11-1552 

8 

12-6046 

12-6014 

9 

14-0503 

14-0475 

10 

15-4961 

15-4935 

11 

16-9419 

16-9396 

12 

18-3877 

18-3855 
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and  the  various  constants  are  defined  above.  The  exact  value  of  the  frequency  parameter 
for  the  first  mode  is  easily  obtained  (for  the  purpose  of  comparison)  by  the  method  of 
reference  [5]  in  which  v  is  expanded  in  a  simple  power  series  or  equivalently,  a  series  of  all 
Legendre  polynomials.  The  exact  values,  to  six  places,  are  obtained  with  only  a  few  terms 
in  the  series  (ten  terms  are  used  to  report  the  results  here).  It  is  evident  from  Table  1  that 
the  simple  expansion,  equation  (22),  yields  results  that  are  quite  accurate  (within  0-3%  for 
all  a).  Thus,  equation  (22)  is  suitable  for  extending  the  expansion  for  the  fundamental 
frequency  in  reference  [4],  which  applies  only  to  rotating  beams  clamped  at  the  axis  of 
rotation  (a  =  0),  to  the  more  general  case  involving  off-clamping  (a  #0). 


6.  EXTENSION  OF  PETERS’  FORMULA 

In  reference  [4],  Peters  obtained  a  formula  for  /u  when  a  =  0  that  was  approximately 
valid  for  all  77.  In  this  section  a  generalization  is  obtained  that  is  approximately  valid  for  all 
rj  and  a. 

For  small  rj  the  formula  will  have  the  form 

M  =/xo(a)  +  Mi(a)h'/2.  (24) 

where  the  expression  for  /z0(a )  is  given  by  equation  (22).  Results  in  reference  [4]  for  a  =  0 
and  in  reference  [6]  for  a  =  0(r?12)  both  suggest  /i] (a  )  =  3/V2.  One  can  therefore  assume 
that  for  all  a 

+  +  (25) 


The  ultimate  justification  for  such  an  assumption  is  that  it  works,  although  it  is  definitely  a 
possible  source  of  error  for  large  a  and  small  77 1/2. 

Before  proceeding  further  one  must  obtain  a  large  77  approximation  for  /x.  For  large  77 
and  a  =  0(1),  let  e  =  I/77  and  y  - n/v-  Thus,  equation  (1)  becomes 

H""-(e/2)[(l  -x2)n']'-ae[(l  - x)u']' -  yu  =  0.  (26) 

For  u  =  «()  +  u,e  +  ■  ■  ■  and  y  =  y<>+  Yif  +  •  ■  •  one  obtains 

e°:uo~  Yo«o  =  0,  (27a) 

f «r-yoMi  =[|(1  - xz)u'o]'  +  a[(l  - x)uo]'  +  yiUa.  (27b) 

The  boundary  conditions  in  equation  (3)  yield  the  first  mode  shape  for  u  : 


u0  =  cosh  (i\x  -  cos  p\x  —  fi  (sinh  0i*-sin  f}\x), 

yu  =  0u  0i  =  1-875104069,  f,  =  0-7340955138.  (28) 

To  suppress  secular  terms,  u0  must  be  orthogonal  to  the  right-hand  side  of  equation  (27b) 
which  yields 


(1  -x2)u! i2  d*  +a 


u 0  d x  =  A  +  Bat, 


(29) 


where 


A  =  1-193336374  and  B  =  1-570878190.  (30) 

The  integrals  in  equation  (29)  are  calculated  as  outlined  in  the  appendix  of  reference  [7], 
Therefore,  for  large  77. 


ix  =  0jT7  +A  +Ba. 


(31) 
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Although  0=0(1)  was  assumed,  it  is  observed  from  reference  [6]  that  when  the  frequency 
parameter  vanishes  for  large  tj  and  a  (negative),  tj  is  linear  with  a.  Hence,  equation  (31) 
may  be  approximately  valid  for  large  77  regardless  of  the  magnitude  of  a.  One  can  now 
assume  a  composite  expansion  for  all  rj  similar  to  that  in  reference  [4]  of  the  form 


/x  =  04\tj  +  *'0 a  +i>i  +  g(a)  +  (3'J2/ir)T}'/2  tan  1  [?r(a  +  ba)/( 3n/2tj’/2)]-  (32) 

For  small  tjI2>  equation  (32)  reduces  to 

ix  st'ocr  +  r'i  +g(a)  +  (3/J2)TjU2  +  O(tj),  (33) 

which  agrees  with  equation  (25).  For  large  17,  equation  (32)  becomes 

H=/3*t)  4- p„a  +  i>i+g(a)  +  a+ba+0(  l/rj),  134) 

so  that,  from  equation  (31) 

A  +  Ba  =  I'd**  +  v\  +  g(a )  +  a  +  bat,  (35) 


a  =  A  -  p\  ~g{a),  b  =  B  -  va.  (36) 

The  final  composite  expansion  is  then 

ix  =p4ri  +  vaa  +  v^+ g(a)  +  (3'^2./ n)ri'12  tan  ’{ n/(3'fl  tj  1/2)[A  -  p,  ~g(a)  +  (B  -  vola]}, 

(37) 
or 

y  =  P*  +  [foa  +  pi  +  g(a)]{}2  +  (3^2/ n)n  tan’1  {ir/2/(3V2)[A  -  Pi  ~g(a)  +  (B  -  p0)q ]}, 

(38) 

with 

P4i  =  12-36236337,  „n=  1-445796491,  e,  =  1-038163742, 

g(a)  =  e’“[A0-  P\  +  (Ao-^i  +  Ai  -pn)a], 

A„=l,  A,  =3/2,  A  =  1-193336374,  5  =  1-570878190, 

y  =  tx/r1  =  (mL4fE[)(w2  +  n2  sin2  0),  Q  =  rf 1/2  =  O'! mL* / El.  (39) 


Table  2 


Frequency  versus  angular  speed  for  a  =  0 


y'12 

equation  (38) 

y'n 

exact  from  fl,  2] 

0 

3-51602 

3-51602 

1 

3-68163 

3-68165 

2 

413710 

4-13732 

3 

479644 

4-79728 

4 

5  58316 

5-58500 

5 

6-44650 

6-44954 

6 

7-35614 

7-36037 

7 

8-29440 

8-29964 

8 

9-25085 

9-25684 

9 

10-2192 

10-2257 

10 

11-1956 

11-2023 

11 

12-1775 

12-1843 

12 

13-1634 

13-1702 
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The  exact  solution  of  equation  (1)  has  been  calculated  in  reference  [2]  by  using  the 
method  of  Frobenius.  It  is  far  more  convenient  for  present  purposes,  however,  to  obtain 
the  exact  solution  from  the  analysis  of  reference  [1]  with  enough  terms  in  the  analysis  to 
ensure  convergence.  Note  that  in  the  finite  element  analysis  of  reference  [1]  the  element 
displacement  is  expanded  in  a  simple  power  series  which  is,  of  course,  equivalent  to  a  series 
of  Legendre  polynomials.  Convergence  to  the  exact  solution  will  only  rapidly  occur  if, 
when  N  +  1  terms  are  taken,  all  terms  of  degree  N  or  less  are  included.  This  explains  the 
relatively  slow  convergence  in  reference  [8]  where  only  odd  Legendre  polynomials  were 
taken.  In  reference  f  1],  only  a  few  terms  are  required  to  obtain  the  exact  solution  to  six 
places  for  the  first  mode  with  one  element.  The  exact  solution  for  the  dimensionless 
frequency  parameter  y 1  2  is  given  in  Tables  2-5  for  dimensionless  rotor  speed  values 
i?  =  0, 1, ....  12  and  dimensionless  off-clamping  parameter  values  a  =0, 01, 1,  and  10. 

Results  from  equation  (38)  are  presented  in  the  tables  for  comparison  with  the  exact 

Table  3 

Frequency  versus  angular  speed  for  a  =  0- 1 


equation  (38) 

3-51602 

3- 70288 

4- 21208 

4- 94054 

5- 80129 

6- 73944 

7- 72347 

8- 73548 

9- 76513 

10- 8063 

11- 8552 

12- 9094 

13- 9675 


exact  from  [1] 


3-51602 

3- 70290 

4- 21225 

4- 94115 

5- 80256 

6- 74142 

7- 72603 

8- 73839 

9- 76815 

10- 8092 

11- 8578 

12- 9116 

13- 9692 


Accession 

m  13  n*M- 

rnc  tV 

Ur-nnno  ::.- 
Ju  nf.il'  ic.it 

By - 

Dlstribut  - . 
!  Aval Ink  :  ■ 


Table  4 

Frequency  versus  angular  speed  for  a  =  lj 


n 

y1'2, 

equation  (38) 

y1'2. 

exact  from  [1,2] 

0 

3-51602 

3-51602 

i 

3-88874 

3-88882 

2 

4-83279 

4-83369 

3 

6-07891 

6-08175 

4 

7-46955 

7-47505 

5 

8-93210 

8-94036 

6 

10-4331 

10-4439 

7 

11-9563 

11-9691 

8 

13-4931 

13-5074 

9 

15-0388 

15-0541 

10 

16-5905 

16-6064 

1! 

18-1464 

18-1625 

12 

19-7055 

19-7215 

u 


18 


D.  H.  HODGRS 


Table  5 

Frequency  versus  angular  speed  for  a  =  10 


n 

y‘'\ 

equation  (38) 

V'\ 

exact  from  [  1  ] 

0 

3-51602 

3-5 1602 

1 

5-38018 

5-40186 

2 

8-79660 

8-88430 

3 

12-5274 

12-6722 

4 

16-3558 

16-5410 

5 

20-2260 

20-4393 

6 

24-1177 

24-3514 

7 

28-0220 

28-2709 

8 

31-9342 

32-1948 

9 

35-8518 

36-1215 

10 

39-7731 

40-0503 

11 

43-6971 

43-9804 

12 

47-6232 

47-9116  j* 

solution.  The  correlation  is  excellent  for  a  =0,0-1  and  1 ,  the  maximum  error  being  of  the 
order  of  01%.  For  a  =  10  the  correlation  is  not  as  good;  the  errors  reach  1-1%.  Note, 
however,  that  the  error  decreases  as  H  becomes  large  or  tends  toward  zero.  Evidently  the 
formula  is  limited  to  a  =  0(1 )  as  assumed  in  the  large  17  expansion,  but  gives  a  reasonably 
good  estimate  even  for  a  =  10.  For  helicopter  and  wind  turbine  blades,  a  is  of  the  order  of 
0- 1 ;  for  gas  turbine  blades  and  satellite  booms,  a  may  be  1  or  larger. 

It  is  evident  that  equation  (38)  provides  an  accurate  and  rapid  estimate  of  the  frequency 
of  a  rotating  cantilever  without  the  need  for  either  digital  computer  programs  or  inter¬ 
polation  of  charts.  It  would  be  desirable  to  include  tip  mass  and  taper  parameters  in  the 
development,  but  the  feasibility  of  such  an  extension  remains  to  be  seen.  The  work  could 
be  extended  to  higher  modes,  but  only  if  one  were  willing  to  work  through  a  great  deal  of 
algebra. 
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